Pair Distribution Function of One-dimensional "Hard Sphere" Fermi and Bose 

Systems 
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INTRODUCTION 

Recent experimental progress in one-dimensional (ID) 
systems renews theoretical interest in such systems. In 
this paper we study the pair distribution function in ID 
"hard sphere" systems. 

Pair distribution function in ID Tonks gas [Hand 
Lieb-Liniger Bose gas 0, H| have been studied pre- 
viously, including the local correlation, -D(O), asymp- 
totic properties at large r, and the zero temperature 
behavior [ISSSBSH, El, EI EDI- The pair distribu- 
tion function for 3D hard sphere boson gas has been cal- 
culated by the authors 1J| . In this work we address the 
problem of pair distribution function in ID hard sphere 
fermion and boson systems with the diameter of the par- 
ticles a^O. 



The pair distribution function is defined as 

D(r 12 ) = p- 2 (^( ri )^t( r2 )^( r2 )^( ri )) 



(1) 



where p = N/L is the density. D (r) is an important 
physical quantity measurable for many liquid systems. In 
this formula ip(r) is the annihilation operator in r space. 

The pair distribution function D(r 12) is also related to 
the diagonal elements of the reduced density matrix [l5| : 

D(r 12 ) = p- 2 Tr[^(r 2 )^(r 1 )p Jv ^t (ri)t/; t (r2)L (2 ) 

where pn — ^o^O' ano - ^0 is the many body wave func- 
tion. 

The meaning of D(r) is: given a particle A at one 
point, the probability of finding another particle B at 
distance r (counterclockwise) is 



pD(r)dr. 

It is obvious that D(r) — * 1 as r — ► 00. also 
pD(r)dr = N - 1. 



(3) 



(4) 



fermions are 2wk/L, where k = -(N-l)/2 to (N-l)/2. 
Here we consider the case that the number of fermions 
in the system is odd, N = 2n+ 1. The normalized many 
body wave function ^(r\, r-i, ■ ■ ■ , r/v) of the system is of 
the form 



■ det 



A 



c 2ri+ 
n-1 

E 2n+1 
n—2 

E 2n+1 



(5) 



c 2n+ 



where €j = exp(i^-rj) / yL. After some calculations, we 
arrive at a compact form for the wave function, 



V ^ ■ V L l<Kj<N 

(6) 

The pair distribution function D(r) for free fermions can 
be exactly obtained from this many-body wave function 
I 



D(r) = 1 - 



sin 2 (^) 



(7) 



Gap Variables For Free Fermions 

In order to prepare for dealing with the problem, 
we introduce gap variables in the case a = as follows. 
Consider a region of the N body Lx Lx ■ ■ • x L coordinate 
system where, modulo L, 



FERMIONS WITH a = 



ri < r 2 < ■ ■ ■ < r N < ri. 



(8) 



Wave Function 9 For Fermions With a = 

First, we consider N free fermions in a one-dimensional 
cyclic interval of length L. The momenta of these 



We shall designate this region as R\. In this region, we 
introduce gap variables {gi} : 



9i = n+i-n, i = 1,2,3, ...,2n + l. 



(9) 
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Obviously 9i — L. With the gap variables the many 
body wave function \P is given by 

* = Q[/(5i)/(52)---/(52„+i)] 

X [/(.9l + .92)/(ff2 + 93) ■ ■ ■ /(<?2n+l + 

X [/(fifl + .92 + .93) ' ' ' f(92n+l +91+ 52)] 



X [f(gi H h fif„) • • • /(S2n+1 + ffl H 1" Sfn-l)] 

(10) 

where /(<7) = sin(7n?/L), and Q is the normalization fac- 
tor given by 



Q 



L : (J_)^[2 JV ( Ar - 1 )/ 2 ]. 



(11) 



'J has (2n + 1) X n factors in total. 
The probability distributions in R\ is 



M 2 dgidg 2 ■ ■ ■ dg N 5(^ .9* - L ) L 



(12) 



with all gt > 0. Now R\ is only one of the (N — 1)! 
regions of the full coordinate space. In each of these 
regions we have the same gap distribution as (12). Thus 
gap distribution probability dP is 



N 



dP=(N- iy.L\y\ 2 d gi dg 2 ■ ■ ■ dg N S(J2 9i ~ L )- (13) 




A 

FIG. 1: Cycle of length L. 

Functions Fi(r) For Free Fermions 

We now evaluate the function D{r) of expression (3) 
using (13). Going from particle A to B, counterclockwise 



in the cycle of length L, as shown in Fig. [1] there may 
be i other particles, with i = 0, 1, • • • , (N — 2). Thus 

pD(r) = P [F {r) + F^r) + ■ ■ ■ + F N _ 2 (r)\, (14) 

where 



„ i+i 

pFi{r)dr = / 5(^2 gj - r)drdP 
J j=i 



(15) 



Thus 



„ i+1 
pF{r) = (N-l)\L / |*| 2 *(j>-r) 

j'=i 



JY 



S ( ^2 g i ~ ( L ~ r )) d 9idg 2 ■ ■ ■ dg N (16) 



where all gj > 0. 

Integrating over dr we get, from (15), 



J pF t {r)dr = J dP = 1. 



(17) 



Thus by (14),/ pD(r)dr = A/ — 1, confirming (4). Outside 
of the interval (0, L) we define F(r) by 



F(r) = for r < and r > L. 



(18) 



Besides (17), i r j(r) has also the following properties: 

(i) It is analytic except at r = and r = L, where F(r) 
and its first derivative are both zero. 

(ii) F i (r)=F N - 2 -i(L-r). 

It is obvious from (16) that NFi(r) is a function of N 
and r/L. In Fig. [2] we plot this function vs r/L for the 
case N = 5. 
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FIG. 2: NFt(r) against r/L for N = 5 and i = 0, 1, 2, 3. 
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FIG. 3: 4 possible configurations which contribute to the pair 
distribution function D(r) for N — 5. The ring represents the 
system and each block stands for a particle with diameter a. 
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F 2 [(r-3a)g 
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FIG. 4: Regions where Fi(r) > for N = 5. Open circles 
are at r — a, 2a, 3a, 4a. Closed circles are at r = L — a, L — 
2a, L — 3a and L — 4a. These 2(N — 1) = 8 circles are where 
D(r) is singular. Each of the four regions has length L — Na. 



FERMIONS WITH a > 

For fermions with a > 0, the full coordinate space is 
again divided into (N — 1)! regions, one of which, R\, is 
defined by the cyclic condition (8) modulo L. We intro- 
duce gap variables in R\ by equations similar to (9): 



9i 
and 



n+i 

5jv = 



- n 
n - 



-a, i = 1, 2, 3, ■ 
tat + L — a 



Obviously 



L- Na. 



,(N-1) (19) 
(19a) 



(20) 



The key point is that "J in terms of the gap variables 
{g} is still given by (10), but with 



1 



1 



Equation (12) now becomes 



jvpjv(jy-i)/2j 



(22) 



M 2 d gi dg 2 ■ ■ ■ dg N S(J2 9* - ( L - Na))L, 



(23) 



and (13) beomes 
dP= (N-iy.L\^\ 2 dg 1 ---dg N 6(Y,9 l -(L-Na)). (24) 



i=l 

The function D(r) is again, as in (14), a sum of (N — 
\)F functions. To illustrate the reasoning we turn to 
Fig. [3] and Fig. 0] for the case N = 5. We have 

D{r) - F [(r-a)e] + F 1 [(r-2o)g] 

+ F 2 [(r-3a)$+F 3 {(r-4a)Z}, (25) 

where ^ = L/(L — Na), and Fi(r) is defined by (16) and 
(18). 

Since Fi(r) is nonzero only in the open interval (0, L), 
for N = 5, 

Fo[(r — a)^] is nonzero only for a < r < L — 4a, 
Fi[(r — 2a)£] is nonzero only for 2a < r < L — 3a, 
F2[(r — 3a)£] is nonzero only for 3a < r < L — 2a, 
Fs[(r — 4a)£] is nonzero only for 4a < r < L — a. 
In Fig. 2]we indicate the regions where Fi[(r — (i+ l)a)£] 
is nonzero for i = to 3. 

We present in Fig. O and Fig. [5] the function D(r) for 
N = 5, and for several values of L/a. 

Returning now to the general case, we conclude 

In the interval (0, L), D(r) is analytic everywhere 
except at (JV — 1) open circles at r = a, 2a, • • • , (N — 
l)a, and at (N — 1) full circles at r = L — a, L — 
2a, • • • , L— (N— l)a. At these singular points D(r) is 
continuous and has a continuous first derivative. 



BOSONS WITH a > 

For the one-dimensional hard sphere boson system, Gi- 
rardeau 0, 0| has shown a Bose- Fermi mapping theorem 
which maps the hardcore boson system to a spinless hard- 
core fermion system, for a > 0. 

*BB(^i,a;2,--- ,xn) = |*fd(2;i,:e2, • • • ,x N )\. 

The pair distribution function is related to the square of 
the wave function. Thus for a > 



fig) = sm(7rg/(L - Na)) 



(21) 



D BE (r) = D FU (r)- 



(26) 
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FIG. 5: (color online) The pair distribution function D(r) as 
a function of distance r for N — 5. The iV — 1 = 4 open circles 
and the N — 1 = 4 closed circles are where -D(Y) is singular. 

For case a = L/(5+), ND(r) is four 5 functions (not 
shown) . 
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FIG. 6: (color online) The pair distribution function D(r) as 
a function of distance r for N = 5. 



